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Abstract 

Predictions for some experimental physical observables in nucleus-nucleus collisions at RHIC 
energies are presented. I utilize the previous suggestion that the retarded, by time about 2- 
3 fm/c, jet absorption in opaque core is a natural explanation of many experimental data. This 
assumption is applied only to the particles with high transverse momentum above 4 GeV/c resulting 
from parton fragmentation. I calculate nuclear modification factor Raa, Raa in- and out of the 
reaction plane, azimuthal anisotropy parameter V2, jet suppression Iaa for the away side jet and its 
dependence versus the reaction plane orientation. The systems under consideration are Au+Au, 
Cu+Cu, at 200 GeV and 62 GeV. Most of numbers are predictions prior to the QM2006. 
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At QM2006 conference we expect that many new experimental data will be shown. I 
already get few requests from different groups to calculate some experimentally observed 



parameters. To be fair,_ 
In previous papers 



resent all these numbers prior to the conference. 
J I propose a simple model to explain experimental data on the 
angular dependence of the nuclear modification factor Raa at high transverse momentum in 
the reaction plane. The model has one free parameter L~2 fm to describe the the thickness 
of the corona area with no absorption and was adjusted to fit the experimental data on 
AuAu collisions at centrality 50-60%. The model nicely describes the Raa dependence for 
all centrality classes. I extract the second Fourier component amplitude, t>2, for high pt 
particle azimuthal distribution and found t>2 should be at the level of 11-12% purely from 
the geometry of the collision with particle absorption in the core. At that time I made a 
prediction for Raa in Cu+Cu collisions at 200 GeV which, as later was found, is in very 
good agreement with experimental data. Physical interpretation of the parameter L could 
be that it is actually retarded jet absorption caused by the plasma formation time T = L/c 
~ 2 fm/c, or at least non-trivial response of strongly interacting plasma to fast moving color 
charge. 

This value of formation time is true for the dilute corona region where density is small. For 
more central interaction region it should be smaller, probably changing like mean distance 
between interaction points or 1/ J\p(x, y)), where p(x,y) is two-dimentional (projected into 
x — y plane) density of Nonary or N part . Extrapolation of T in a such way gives values of 
the order of 0.7-1 fm/c in the center of collision zone. 

As usual, nuclear modification factor Raa is defined as: 



(1/JVeg) d 2 N A+A /dp T dr] 
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where (Nbi nar y) is a number of binary nucleon-nucleon collisions at particular centrality class. 

I use a Monte Carlo simulation of nucleus-nucleus collisions based on the Glauber ap- 
proach with Woods-Saxon density distribution. Assume that jets are produced in accordance 
with the Neon distribution in the x — y transverse plane. I neglect longitudinal expansion 
for the first 2-3 fm/c. Let all jets fly in all possible directions (some of them will move from 
the interaction zone determined by the envelop of Woods-Saxon radii, some will go into this 
zone) for some time T. If after that time jet will be still inside the interaction zone, it will 
be completely absorbed. If not, jet leaves the zone without interaction. We found from ex- 



perimental data that this time should be about 2.3 fm/c at 200 GeV and 3.5 fm/c at 62 GeV 
center of mass beam energy. Absorption can not happen instantly (or in the infinitely thin 
layer), plus there should be statistical fluctuations. To take these into account I smooth the 
cut edge: an arbitrary weight function in the form of a Fermi distribution with diffuseness 
parameter a was applied: 

Wei9ht{l) = l + expUT-iyaY (2) 
The value of this function changes from 0.1 to 0.9 at Al=4.4*a around T. Parameter a 

varied from 0.01 to 0.5 fm/c. It was found that the results for Raa and v% do not change 

much for this parameter range. By default it was chosen to be a=0.2 fm/c. The results of 

calculation for Raa are presented in Fig. CD All values for Raa and v% are presented in tables 

at the bottom of this preprint. It is worth to mention, values of Raa and Vi for mesons and 

baryons at high p? should be the same. Obviously, Raa is flat versus pt if no other physics 

process is "switched on" at high p T , except retarded absorption in opaque core considered 

here. 

For di-jets, another parameter, Iaa, was introduced in Ref. [51]. For away side jet Iaa 
is the ratio of particle yield per trigger in the particular momentum range of trigger and 
associate particles, to the similar yield in p+p collisions. To calculate this value I perform the 
same procedure like for Raa, except that in addition the location of away side jet after time 
T should be taken into account. Away side jet may not go exactly in opposite relative to the 
trigger jet direction, but is distributed around 180 degree angle. It introduces another critical 
parameter for the calculation, which should be taken from the experimentally observed width 
of the away side jets in azimuthal angle. I use simple Gausse distribution around 180 degrees 
for the away side jet. Depending on the momentum range of the associate particle, width 
of the away side jet varies from 0.5 radians for the low momentum range to 0.2 radians at 
momentum above 5-6 GeV/c. 

Results of calculation for Iaa in- and out of plane for Au+Au collisions at 200 GeV are 
presented in Fig. [2j One can see the main feature of the surface biased scenario: in the 
most central events there are more di-jets out of the reaction plane than di-jets in the plane. 
This is opposite to the Raa reaction plane dependence or to the scenario which assumes 
only "punch-through" di-jets. In mid-central events the situation is mixed, and in peripheral 
collisions there are more in-plane di-jets. Wide centrality bin of 20 to 60%, commonly used 



in experiments, corresponds to N part =125 if number of di-jets is weighted by the number of 
trigger particles. Resuts for Cu+Cu collisions are shown in Fig. [31 

I want to emphasize that the role of event by event fluctuations should be treated by a 
better way than I do it here by parameter a. It is known that nucleon position fluctuations 
in the colliding Cu+Cu nuclei can significantly increase the observed v% in the most central 
collisions [6]. This is why static models, including our model presented here, underestimate 
t>2 in the central Cu+Cu collisions. The same should be true when we consider di-jets in 
Au+Au or Cu+Cu: local density fluctuations induce positive correlation between surviving 
probability of the trigger jet and the away side jet. Again, such fluctuations could increase 
the observed value of Iaa-, especially in the central collisions. 

Few more comments: as I pointed out in the previous paper [2|, if we see a significant 
suppression of non-photonic electrons p|, one can expect values of t>2 for such electrons at 
high pt to be equal to V2 for high p? hadrons. Another point: suppression of the away 
side jets in direct gamma-hadron correlations, Iaa, should be equal to the Raa of inclusive 
hadrons at about gamma energy. Here I ignore the difference in quark and gluon fragmenta- 
tion functions. The reason, that this may change the result, is that hadrons in the away side 
jets for the direct gamma-hadron correlations are mostly from quark fragmentation, while 
inclusive hadrons below 10 Gev/c are dominantly from gluon fragmentation. 

All numbers for Au+Au collisions are presented at 200 GeV in Table (TJ at 62 GeV in 
Table [TT1 Numbers for Cu+Cu collisions at 200 GeV are in Table [TTTl and at 62 GeV are in 
Table ED 
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Figures 




FIG. 1: Calculated Raa for Au+Au and Cu+Cu collisions at 200 GeV and 62 GeV versus the 
number of participant nucleons, N part . Solid curves are for 200 GeV, dashed are for 62 GeV. The 
circles are experimental it data for 200 GeV Au+Au collisions integrated for j>t>4 GeV/c j^j. 
The triangles are data for 200 GeV Cu+Cu collisions of pr>7 GeV/c Only statistical errors 
are shown. 
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FIG. 2: Calculated Iaa in- and out of plane for Au+Au collisions at 200 GeV versus the number 
of participant nucleons, N par t. Width of the away side jet is <r=0.22 radians. Small triangles are 
estimation for in plane, squeres- for out of plane. The circles with error bars are data from j|. 
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FIG. 3: Calculated Iaa for Cu+Cu collisions at 200 GeV versus the number of participant nucleons, 
Npart- Width of the away side jet is a=0.22 radians. Circles are estimation for in plane, triangles 
- for out of plane. 
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TABLE I: Results of calculations for Au+Au at 200 GeV. In all calculations a=0.2 fm/c and 
o"=0.22 radians. 



cent, % N p art 


■f^ R-in Rout 


V2, % 


< Iaa > 


In lout 


0-5 


353 


0.195 0.213 0.177 


4.71 


0.145 


0.0882 0.202 


5-10 


300 


0.235 0.281 0.211 


7.06 


0.198 


0.135 0.262 


10-15 


255 


0.29 0.348 0.246 


8.64 


0.259 


0.204 0.314 



15-20 215 0.335 0.417 0.282 9.68 0.34 0.317 0.364 
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20-25 182 0.39 0.483 0.316 10.47 0.448 0.488 0.407 



25-30 152 0.45 0.563 0.362 10.87 0.575 0.687 0.463 

30-35 126 0.505 0.632 0.40 11.23 0.659 0.809 0.508 

35-40 103 0.565 0.705 0.447 11.19 0.727 0.897 0.557 

40-45 84 0.635 0.779 0.512 10.37 0.783 0.947 0.62 

45-50 67 0.715 0.848 0.589 9.03 0.833 0.974 0.69 

50-55 52 0.79 0.902 0.681 7.0 0.877 0.99 0.763 

55-60 40 0.84 0.958 0.847 5.53 

cent, % Np ar t R Ri n Rout V2,%< Iaa > hn lout 

TABLE II: Results of calculations for Au+Au at 62 GeV. 

Cent, % Npart R Rin Rout «2, % 

0-5 345 0.333 0.364 0.302 4.65 

5-10 296 0.402 0.461 0.344 7.29 

10-15 250 0.468 0.553 0.383 9.06 

15-20 211 0.532 0.639 0.425 10. 

20-25 178 0.591 0.719 0.462 10.9 

25-30 149 0.656 0.801 0.512 11. 

30-35 123.8 0.707 0.861 0.553 10.9 

35-40 102 0.756 0.911 0.601 10.3 

40-45 82.9 0.805 0.948 0.663 8.83 

45-50 66.2 0.851 0.971 0.731 7.05 

50-55 51.9 0.895 0.984 0.805 5.01 

55-60 40 0.922 0.99 0.854 3.69 

Cent, % Np art R Rin Rout V 2 , % 

TABLE III: Results of calculations for Cu+Cu at 200 GeV. 
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cent, % Npart R R in R out v 2 , % < Iaa > hn hut 



0-5 110 0.37 0.40 0.34 3.4 0.194 0.136 0.251 



5-10 


95 


0.435 0.48 0.39 


5.1 


0.263 


0.203 0.324 


10-15 


82 


0.51 0.57 0.45 


6.0 


0.346 


0.31 0.382 


15-20 


70.5 


0.585 0.66 0.51 


6.3 


0.482 


0.497 0.467 


20-25 


60 


0.655 0.74 0.57 


6.4 


0.581 


0.637 0.526 


25-30 


51 


0.7 0.79 0.61 


6.4 


0.645 


0.723 0.568 


30-35 


43 


0.775 0.86 0.69 


5.4 


0.733 


0.817 0.648 


35-40 


36 


0.825 0.9 0.75 


4.3 


0.784 


0.865 0.704 


40-45 


29 


0.87 0.93 0.81 


3.3 


0.833 


0.903 0.762 


45-50 


24 


0.91 0.95 0.87 


2.1 






50-55 


19 


0.95 0.97 0.93 


0.9 






55-60 


15 


0.98 0.98 0.98 










Cent, % N part R Rin Rout V 2 , % < IaA > lin lout 



TABLE IV: Results of calculations for Cu+Cu at 62 GeV. 



cent, % N part 


-Rin ^out 


V2, % 


0-5 


104 


0.598 0.638 0.559 


3.28 


5-10 


92 


0.672 0.737 0.607 


4.82 


10-15 


80 


0.736 0.817 0.656 


5.48 


15-20 


68 


0.8 0.888 0.712 


5.48 


20-25 


57 


0.845 0.931 0.759 


5.08 


25-30 


48.6 


0.874 0.954 0.793 


4.61 


30-35 


40.6 


0.91 0.975 0.846 


3.52 


35-40 


34 


0.932 0.982 0.882 


2.7 


40-45 


28 


0.952 0.989 0.915 


1.93 



Cent, % Npart R Rin Rout V 2 , % 
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